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ABSTRACT 

The stability of a quasicrystalline structure, recently obtained in a molecular-dynamics simulation 
of rapid cooling of a binary melt, is analyzed for binary hard-sphere mixtures within a density- 
functional approach. It is found that this quasicrystal is metastable relative to crystalline and fluid 
phases for diameter ratios above 0.83. Such trend is partially reversed for lower diameter ratios, 
since the quasicrystal becomes stable with respect to the crystal but does not reach a coexistence 
with the fluid. 
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§1. INTRODUCTION 



It is well known that quasicrystalline phases are obtained from alloys of at least two metallic el- 
ements, mostly by processes of rapid solidification of the melt. However, the question of how 
quasicrystals form is far to be understood from the theoretical viewpoint (Steinhardt and Di Vin- 
cenzo 1991, Janot 1994). In this respect, an important step towards this goal could be given by 
a systematical study of numerical simulations at which a quasicrystalline phase is formed from a 
binary fluid. Unfortunately, this program is difficult to carry out since simulations of such systems 
are very expensive due to the large number of atoms which are necessary to display the quasiperi- 
odic features. Recently, however, it was reported a first effort of a molecular dynamics simulation 
of quasicrystal solidification from the melt of a diatomic three-dimensional system (Roth, Schilling 
and Trebin 1995). In fact, it was shown that a binary Lennard- Jones fluid, at adequate cooling 
rates, indeed solidifies into a phase of icosahedral long-range order close to a quasicrystalline struc- 
ture known as the truncated icosahedral binary model (TIBM). This is assumed to be a simplified 
version of the Henley-Elser model for the icosahedral structure of (Al,Zn)4gMg32 (Henley and Elser 
1986). A previous investigation by means of relaxation simulations had shown that the TIBM was 
at least strongly metastable for T=0, a feature which was mainly interpreted due to the fact that 
such structure comprises a very dense packing of spheres of two sizes (Roth, Schilling and Trebin 
1990). 

In the present paper we shall investigate the stability of a TIBM composed by hard-spheres 
by means of a density- functional theory (Lutsko and Baus 1990), which was previously applied to 
a one-component hard-sphere (HS) quasicrystal (Cataldo and Tejero 1995). HS solids may serve 
as suitable reference systems for mean-field- or perturbation-theory treatments of more realistic 
potentials (Daanoun, Tejero and Baus 1994, Achrayah and Baus 1998). Recently, Denton and 
Hafner ( 1997a, b) have studied a model of monatomic quasicrystal interacting via effective metallic 
pair potentials, by means of a combination of classical density-functional theory and thermodynamic 
perturbation theory, which takes as its reference system the above mentioned one-component HS 
quasicrystal. 

Closely related to freezing and thermodynamic stability properties of HS solids is the problem 
of packing. This is a well studied subject for identical hard spheres, but less is known about packing 
properties of mixtures of spheres or disks of two different sizes, see e. g., Likos and Henley (1993). 

The possibility of quasicrystalline metastable states for one- and two-component HS solids has 
been recently explored within density-functional theory in the range of diameter ratios above 0.85 
(M c Carley and Ashcroft 1994). Then it was found that such phases have higher free energies 
than the crystalline or liquid states. We shall see that a similar conclusion stands for the TIBM 
for diameter ratios above 0.83. However, this trend is partially reversed for lower diameter ratios, 
since the quasicrystal becomes stable with respect to crystal but remains metastable with respect 
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to fluid. 

In the following Section we shall describe the TIBM sphere packing, showing that a small shift in 
the original position of large spheres makes the maximum packing fraction to become competitive 
with respect to HS crystal structures of similar composition. In Section 3 we shall give a brief 
account of the density-functional approach to the calculation of quasicrystal free energy. Finally, 
our results and conclusions will be presented in Sections 4 and 5, respectively. 

§2. TIBM SPHERE PACKINGS 

Let us consider a decoration of the 3D Penrose tiling at which all vertices and middle edges of 
both rombohedral cells are decorated by identical spheres. Thus, if we take rombohedral edges as 
unity, the maximum diameter turns out to be 0.5. On the other hand, each prolate rombohedron 
can accomodate two larger spheres located symmetrically at both sides of the center of its long 
diagonal (see fig. 1). If the location of such spheres divides the diagonal into three parts of ratios 
r~ 2 : r~ 3 : r~ 2 , with r = (v5+l)/2, the maximum diameter (given by the separation of the centers 
of both spheres) is 0.563, yielding a diameter ratio 0.888 and a maximum packing fraction of 0.580. 
Here it is convenient to define a notation for the HS parameters namely, a = o~i/o~2, X2 = N2/N, 
and 77 = ^[(1 — X2)a 3 + X2]yO~2, denote respectively, the diameter ratio, the concentration of large 
spheres, and the packing fraction. 

Now, taking into account that in the Henley-Elser model the above ideal position of the larger 
spheres was chosen essentially for simplicity (Henley and Elser 1986), we may allow a greater 
separation between both spheres in order to increase the maximum diameter. In fact, by this 
procedure the maximum packing fraction can be improved to r\ max = 0.648, for a diameter ratio 
a = 0.798 when large spheres of adjacent rombohedra touch each other (in addition to the contact 
between both large spheres of each rombohedron). 

In order to measure the relevance of the above packing fractions, we may consider a binary 
crystal of similar composition namely, a fee lattice plus a decoration of small spheres on faces 
and large spheres on vertices, which has a concentration X2 = 0.25 close to the TIBM value 
(0.236). In addition, we may compare to a disordered- fee structure in which both kind of spheres 
are distributed irregularly over the sites of a fee lattice. This disordered crystal allows us to select 
the exact composition of the TIBM. 

In table |l] we give the packing fractions of the above structures for three selected values of the 
diameter ratio. We note that the diameter ratio 0.798 yields the better quasicrystal performance 
with a packing fraction slightly higher than the crystal one (the packing fraction of the disordered 
crystal is only relevant at the higher diameter ratio). 
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Table 1: Maximum packing fraction of the TIBM quasicrystal (q), the crystal (c), and the disordered 
crystal (d), for three selected values of the diameter ratio a. 



a q c d 

0.888 0.580 0.683 0.572 

0.833 0.618 0.658 0.502 

0.798 0.648 0.643 0.462 

§3. CALCULATION OF THE HELMHOLTZ FREE ENERGY 

The above calculation of packing fractions may be regarded as an "order-zero" of a stability analysis 
of the competing HS structures. In order to get a much more definite analysis, in what follows we 
shall resort to a density-functional approach, the so-called generalized effective liquid approximation 
(Lutsko and Baus 1990), which has proven to provide a very accurate description of one-component 
HS solids. In particular, the application of this formalism to a monatomic HS quasicrystal of high 
packing fraction (Henley 1986), was recently reported (Cataldo and Tejero 1995), showing that it 
is thermodynamically unstable (i. e., it does not have a coexistence with the fluid), in agreement 
with previous studies (M c Carley and Ashcroft 1994). Thus, our present treatment will be based 
upon a straightforward generalization to mixtures of the above approach. In fact, following Xu and 
Baus (1992), we write the Helmholtz free energy per particle of a HS binary solid as the following 
functional of both partial local number densities Pi(r) (i = 1, 2), 

f[pi,P2] = -^E/ dr«WPn(A?Pi(r)) - 1] + f ex [pi, P2}, (1) 
i=i 

where ks denotes Boltzmann's constant and Aj the thermal de Broglie wavelenght of species i. The 
excess free energy per particle of the solid f ex is approximated by the one of an effective HS binary 
fluid of the same concentration, fex/^BT — ^(t/), whose packing fraction 17 is obtained from the 
solution of a system of two coupled differential equations (Tejero and Cuesta 1993), 

_ z(\) - 4>(fj(\)) 

z'(X) = d>(r/(A)) (3) 

as 77 = r)(X = 1). In the above equations, the prime denotes function derivative, the initial conditions 
are 7/(0) = z(0) = 0, and 
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HvW) = ~^ [drf drXEft( r )^( r/ )C^(|r-r'|,7/(A)), 
J J i=i 7=1 



(4) 
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where Cfj denotes the Percus-Yevick approximate direct correlation function of a binary HS fluid 
mixture (Lebowitz 1964). The HS density profiles are parametrized in terms of Gaussians centered 
at the sites {RW} occupied by species i: 

^(r)=(-) 3/2 E^ (r - R<!))2 * = 1 > 2 - (5) 
W 5(0 

Replacing (||) in (||), and performing the angular integrations, we obtain^] 

, 2 2 N t n 3 



being 



S(r, aij,nj) 



*=1 J=l R(0 R(i) 
1/2 



exp(-ay(r - ry) /2) - exp(-ay(r + ry) /2) , (7) 



where Ojj = 2ajOj/(aj+a :) ') and = |R^ — Rw)|. The calculation of (||) deserves some discussion. 

First, we note that the integral in r can be analytically computed and, in the case of periodic 

structures, the double lattice sum reduces to a simple one running over spherical shells of sites 

centered around an arbitrary site. On the other hand, the lack of periodicity of the quasicrystal 

does not allow this simplification, making necessary to resort to another method. In a previous 

work (Cataldo and Tejero 1995), it was proposed a method of evaluation of quasilattice sums that 

can be applied to the present case. We refer the reader for details to such work, but in a brief 

summary the method consists in a generalization of the above mentioned shell summations of 

periodic lattices to a situation in which each site produces a different "shell" sum due to variations 

in the surroundings, as is the case of a quasilattice. So, the double sum consisting in the average 

of such "shell" sums over an infinite number of sites is approximated by an average over the sites 

within a finite size quasicrystal "sample" . Here it is important to notice that in the case of a crystal, 

this procedure yields the correct value only if the "sample" possesses the exact composition of the 

infinite crystal (note that it is subjected to variations due to boundaries). So, we have adopted this 

exact composition prescription for the quasicrystal "samples" . 

The Gaussian width parameters {01,02} of (||) are determined so as to minimize the total 

free energy of (|l]). Fig. 2 shows the effect of quasilattice sums on the value of quasicrystal free 

energy; we see that it suffices to work with "samples" of N2 > 1000 large spheres in order to reduce 

convergence errors to less than 0.1%. In fact, the computation time needed for evaluating each 

(0:1,0:2) f ree energy value is around 11' CPU in a workstation Alpha DEC 3000 for a "sample" of 

4000 large spheres. Fortunately, the location of free energy minima in the 01,02 plane is rather 

1 In case of a disordered structure, instead of (|H|) we must consider the alternative expression, pi(r) = 
Xj (-^-) 3 ^ 2 y"] p exp[— Qi(r — R.) 2 ] (i = 1,2), where the summation runs over all lattice sites. Accordingly, the 
expression (^|) has to be modified by the replacements, X*r(0 — * Xi I]r- 
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insensitive to the "sample" size, so we could work with small "samples" (N2 ~ 100) in locating 
minima and, at a second stage, a minimum value of free energy was computed for a large "sample" 
(N2 ~ 4000) within a small grid of typically 20 (0:1,0:2) points. 

An alternative treatment of the problem of quasilattice sums consist in replacing the quasicrystal 
by a suitable crystalline approximant (Krajcf and Hafner 1992, Goldman and Kelton 1993). Such 
structures are crystals with large unit cells where the arrangements of sites closely approximate 
the local structure of the quasicrystal. Instead of an aperiodic tiling of space by prolate and oblate 
rombohedra (Penrose tiling), approximants are based upon tilings by the same rombohedral units 
which are periodically repeated in the form of large cubic unit cells. Once such tiling is defined, the 
decoration of both kind of rombohedra with large and small spheres can be done in the same way 
as in the quasicrystal. Lattice sums over the approximant are then computed by taking any of its 
cubic unit cells as a "sample" in the above explained sense. It is clear that approximants having 
larger unit cells are expected to yield free energy values closer to the quasicrystal ones. However, 
for HS potentials only the most local structure in the neighbouring of each particle should be of 
importance, and then one expects that the smallest approximants would be adequate to obtain 
acceptable results. We have found that this is indeed the case; one of the smallest approximants 
is the so-called '1/1' cubic structure (Elser and Henley 1985), which has 20 prolate and 12 oblate 
rombohedra per unit cell. It is slightly richer in large spheres (+0.84%) than the quasicrystal and 
thus has a 0.03% (0.11%) higher packing fraction for o=0.888 (a=0.798). Regarding its free energy, 
it turns out to be always slightly lower than the quasicrystal one, namely, for o=0.888 the free 
energy differences vary from 0.13% (nearly within convergence errors of the quasicrystal values) at 
low packing fractions, up to 0.21% at high packing fractions, while for a=0.798, the values increase 
from 0.21% up to 0.41%. From all these figures, it becomes clear that such differences in free energy 
are most probably due to the slightly higher packing fraction of the 1/1 approximant. 

§4. RESULTS 

In fig. 3 we represent the free energy of the different phases for a = 0.888, note that the quasicrystal 

turns out to be always metastable with respect to the remaining phases. On the other hand, both 

crystals show an exchange of relative stability with respect to the fluidQ, which is consistent with a 

freezing into the disordered solid as expected (M c Carley and Ashcroft 1994, Xu and Baus 1992). 

Another feature of solid phases concerns to the degree of localization of the spheres at the 

corresponding sites. This is measured by the Gaussian widths which are shown in table ||. Note 

the particular behavior of small spheres in the quasicrystal, they present an appreciably lower 

localization than the large or small spheres of the remaining cases. 

2 The fluid phase free energy was computed from the BMCSL equation (Boublik 1970, Mansoori, Carnahan, 
Starling and Leland 1971). 
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Table 2: Gaussian width (ccj) 2 in units of the respective HS diameter Oi (diameter ratio a = 0. 
for the quasicrystal (q), the crystal (c), and the disordered crystal (d). 



V 


1 






(02)7* 


MP 




0.52 


0.20 


0.11 


0.12 


0.11 


0.13 


0.14 


0.53 


0.17 


0.10 


0.11 


0.094 


0.12 


0.12 


0.54 


0.16 


0.093 


0.10 


0.085 


0.11 


0.11 


0.55 


0.15 


0.086 


0.097 


0.078 


0.10 


0.10 


0.56 


0.14 


0.080 


0.090 


0.071 


0.092 


0.094 


0.57 


0.13 


0.075 


0.082 


0.064 


0.085 


0.090 


0.58 


0.12 


0.070 


0.076 


0.058 






0.59 






0.070 


0.053 






0.60 






0.064 


0.048 






0.61 






0.059 


0.044 






0.62 






0.054 


0.039 






0.63 






0.049 


0.036 
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An additional test on quasicrystal stability may be performed by taking the distance D between 
the centers of both large spheres of each rombohedron as a variational parameter. In fact, in fig. 4 
we find free energy minima for D values which are approximately equidistant from the original D in 
the Henley-Elser model and the one that allows to maximize the packing fraction (both distances 
are indicated by dotted vertical lines). It is interesting to note that the increase of free energy at 
both sides of minima is accompanied by a loss of particle (small and large spheres) localization. 
Moreover, for each packing fraction there exists a critical maximum separation (star points in fig. 
4) beyond which the free energy minimum in ai,«2 vanishes, yielding thus an end of mechanical 
stability. In order to understand such behavior, we first note from fig. 4 that the critical D increases 
with rj, and then for a fixed D, say 0.66, the quasicrystal becomes mechanically unstable below 
certain density, rj ~ 0.53 in this case. In other words, the bounds of mechanical stability in fig. 4 
denoted as star points, were to be expected as arising from the familiar fact that any solid becomes 
unstable below certain density. 

Let us now repeat our analysis for the lower diameter ratio a = 0.798. Fig. 5 represents the free 
energy of quasicrystal (for D = 0.627), crystal, and fluid. We see that the improved quasicrystal 
packing fraction has contributed to stabilize such phase with respect to the crystal, but it remains 
unstable as compared to the fluid. Note that the crystal free energy is very close to the quasicrystal 
one at ry = 0.53, but unfortunately it was impossible to perform a reliable calculation of the crystal 
free energy for r] < 0.53 since the ai, 0.1 minimum at those low densities began to suffer an overlap 
with a spurious minimum arising from the flaw of Percus-Yevick approximate direct correlation 
function of the effective fluid at solid-like densities^. On the other hand, we have found that the 
crystal minimum in a\,a2 vanishes for r\ > 0.59, this unexpected behavior seems to be caused by 
some kind of high-density flaw of our density- functional theory^]. Nonetheless, the crystal density 
range shown in fig. 5 is enough to perform our comparison of relative stabilities. The Gaussian 
widths of both solid phases are shown in table 3, again we observe that the small spheres in the 
quasicrystal are less localized than the large ones, but both spheres are more localized than the 
corresponding values of table 2. Regarding the crystal, note the anomalous behavior of the last 
column; the loss of localization of large spheres at increasing densities finally yields a vanishing 
minimum in the 02 coordinate. Finally in fig. 6 we test again the quasicrystal stability against 
variations of the distance D between large spheres. We observe features similar to those discussed 
from fig. 4 above, despite the shorter range of D due to a longer large sphere diameter. 

Concluding our results, we note that the exchange of relative stabilities between crystal and 

quasicrystal (cf. figs. 3 and 5) occurs at a diameter ratio 0.83 (see table 1) when the free energies 

of both solids differ within the order of quasicrystal convergence errors, along the whole range of 

3 Such spurious minima were discussed by Lutsko and Baus (1990) for one-component HS crystals. 
4 Recently it has been shown that the density- functional theory employed by M c Carley and Ashcroft (1994) is 
unable to predict the existence of certain high density, mechanically stable, HS solids (Tejero 1997,1998). 
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Table 3: Same as table 2 for a diameter ratio a = 0.798. The disordered crystal is absent due to 



its low packing fraction (see table 1). 

_i _i _i _i 

V 2 ("2)9 2 («i)c 2 (a 2 ) c 2 

0.52 0.16 0.085 

0.53 0.14 0.078 0.15 0.14 

0.54 0.13 0.072 0.14 0.13 

0.55 0.12 0.067 0.13 0.12 

0.56 0.11 0.061 0.12 0.11 

0.57 0.10 0.057 0.12 0.12 

0.58 0.096 0.053 0.11 0.13 

0.59 0.088 0.049 0.11 0.15 

0.60 0.081 0.046 



0.61 0.075 0.043 

0.62 0.069 0.041 

0.63 0.064 0.039 

0.64 0.059 0.037 

packing fraction. 

§5. SUMMARY AND CONCLUSION 

We have analyzed the stability of a HS binary quasicrystal by means of a density-functional method. 
We have first focused upon a TIBM structure with a diameter ratio 0.888 that maximizes the 
packing fraction. It was found that this HS structure turns out to be thermodynamically unstable 
i. e., it has a higher free energy with respect to crystals and fluid of similar composition, being a 
disordered crystal the solid phase of lowest free energy. Next we have studied the quasicrystal free 
energy dependence upon the distance between nearest large spheres, finding a minimum at a slightly 
more distance than the original TIBM one. In addition, we have found that a small increase of 
such distance allows to improve the maximum packing fraction to a competitive value with respect 
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to a crystal of similar composition. In fact, it was found that for a diameter ratio around 0.8, the 
quasicrystal is stable with respect to such crystal, but it remains unstable as compared to the fluid. 

In conclusion, as was pointed out by M c Carley and Ashcroft (1994), it appears that crystalline 
states are preferred for diameter ratios above 0.85 but, as was shown in the present paper, a small 
decrease of diameter ratio to 0.83 leads to a favored quasicrystalline state which, however, does 
not reach a coexistence with the fluid. Such thermodynamical unstability is most probably due 
to the assumed HS interparticle potential, in contrast to the Lennard-Jones one of the simulation 
performed by Roth et al. (1995). 
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Figure 1: Decoration of a prolate rombohedron containing two large spheres located symmetrically 
at both sides of the center of its long diagonal denoted by a dotted line. The maximum diameter of 
small spheres can reach half a rombohedron edge. The separation of both large spheres is bounded 
by the contact with other large spheres located in adjacent rombohedra, or eventually by the contact 
with the small spheres denoted as 1, 2, 3. 
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Figure 2: Free energy per unit volume = rj[f /ksT — 31n(A2/o"2) + 1], of spherical quasicrystal 
"samples" of different sizes containing N2 large spheres. The HS parameters are a = 0.798, rj = 0.55, 
and the distance between both large spheres of each prolate rombohedron is 0.627. 
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Figure 3: Free energy per unit volume (ft = ry[/ /ksT — 31n(A2/o"2) + 1], vs packing fraction 77, for 
quasicrystal (diamond points), crystal (square points), disordered crystal (circle points), and fluid 
(solid line). 



5.0-1— r 




Figure 4: Quasicrystal free energy per unit volume (ft, vs the distance D between both large spheres 
of each rombohedron, for different packing fractions. The circle points (left end of curves) indicate 
that both spheres are close to touch each other, while the star points (right end of curves) indicate 
an end of mechanical stability. The dotted vertical lines correspond to D = 0.563 and D = 0.627, 
the diameter ratio is a = 0.888. 
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Figure 5: Same as fig. 3. The disordered crystal is absent due to its low packing fraction (see table 
1)- 
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Figure 6: Same as fig. 4 for a diameter ratio a = 0.798. The circle points at the right end of upper 
curves indicate that each large sphere is close to touch the small spheres denoted as 1, 2, 3 in fig. 
1. 
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